Observable core response in neutron star spin glitches 
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The superfluid core of a neutron star is usually assumed to corotate with the crust over timescales 
longer than minutes. I show that the interaction between the neutron superfluid and the type II 
superconductor of the outer core increases the coupling time to weeks or longer. I suggest that 
observed post-glitch response over timescales of weeks to years represents recovery of the outer core. 
Spin glitches could originate in either the inner crust or the outer core. 
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Introduction. The problem of the origin of spin glitches 
in neutron stars (NSs) has remained unsolved since the 
first glitch was observed in the Vela pulsar in 1969. This 
problem is of considerable interest, as an understand- 
ing of the glitch phenomenon would offer insights into 
both the dynamical and ground-state properties of mat- 
ter above nuclear density. To make progress on the glitch 
problem, it is crucial to identify which components of the 
liquid interior corotate with the crust, store and release 
the angular momentum that drives glitches, and produce 
the observed post-glitch response that occurs over days 
to years (see [lj for examples). 

The outer core of a NS was predicted to contain su- 
perfluid neutrons and superconducting protons long ago 
(0; for a recent review, see Q). The neutron flow around 
the vortices that thread the rotating superfluid entrains a 
proton mass current, giving the vortices very strong mag- 
netization [4j . Electrons scatter efficiently with the mag- 
netized vortices, enforcing corotation between the neu- 
tron fluid and the proton-electron fluid over a timescale 
of 400 - 10 4 P, where P is the rotational period in sec- 
onds [5j. For typical rotation periods of < 0.1 s, the cou- 
pling time is minutes. Corotation between the charged 
core fluid and the crust is maintained by a combination 
of Ekman pumping and magnetic stresses The en- 
tire core is thus usually assumed to corotate with the 
crust. If this assumption is correct, glitches can arise 
only in the inner-crust superfluid that coexists with the 
solid lattice, where the pinning of vortices to the nuclear 
lattice produces metastable rotational states 0. Post- 
glitch response is then naturally attributed to relaxation 
of the inner-crust superfluid (e.g., [8|]). Ignoring the natal 
magnetic field of the NS, Ref. [9| argued that magnetic 
flux tubes of the outer core could form clusters around 
vortices, which could increase the coupling time to ob- 
servable timescales of days or longer. 

In this Letter, I revisit the issue of core coupling con- 
sidering an essential effect of type II superconductivity in 
the outer core, that of the pinning of vortices to magnetic 
flux tubes associated with the star's natal magnetic field. 
As the protons condensed when the star was young, the 
very high electrical conductivity of the relativistically- 



degenerate electrons prevented Meissner expulsion of the 
core magnetic field [10( , and the field formed a dense tan- 
gle of magnetic flux tubes with which vortices interact. 
This primarily magnetic interaction pins the vortices to 
the magnetic tangle which is frozen to the superconduct- 
ing fluid. I show that the coupling time between the neu- 
tron and charged fluids is days to years, so that glitch re- 
covery involves the response of the supra-nuclear- density 
matter of the outer core, where pinning occurs through 
a very different mechanism than in the inner crust. The 
problem addressed here was also considered by Ref. [Il[ , 
who concluded that the relaxation time is << 1 d; the 
analysis given here gives a much longer longer timescale. 

Vortex drag. The equation of motion of a straight 
vortex moving at velocity v with respect to the proton- 
electron fluid is given by balance between the Magnus 
force on the vortex and the drag force. In the rest frame 
of the charged component: 

p s K x (v - v s ) - r\v = 0, (1) 

where v s is the flow velocity of the ambient superfluid, p s 
is the superfluid mass density, k = h/2m n is the vortic- 
ity quantum where m n is the neutron mass, k is directed 
along the vortex, and r\ is the drag coefficient. In cylindri- 
cal coordinates (r, 4>, z) , with z and n along the rotation 
axis, the vortex velocity is given by |12| 

v = v s Qsin20 d f + cos 2 9 d ^j , (2) 

where the dissipation angle 0d is given by 

tan6» d = — . (3) 
p s n 

The vortex moves at an angle 9d with respect to v s . En- 
trapment of the neutron and proton mass currents when 
both species are superfluid endows a vortex with a mag- 
netic field of B v ~ 10 14 G Electron scattering with 
the strongly-magnetized vortex cores gives [f| 

tan0 rf = (Qst)- 1 ~ 10" 3 , (4) 

in the outer core, where £l s is the spin rate of the su- 
perfluid and t is the dynamical relaxation time between 
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the neutron fluid and the charged fluid. While the dissi- 
pation angle is much smaller than unity, it is still large 
enough that the neutral and charged fluids are coupled 
over a timescale of minutes. 

Pinning of vortices to flux tubes. The magnetic field 
in a flux tube is -B$ ~ 10 15 G [J|. The magnetic inter- 
action energy between a vortex and a flux tube E p is of 
order B V B$,V, where V is the overlap volume, and has 
been estimated to be ~ 100 MeV by a number of authors 
Til 14 1 . The angle-averaged interaction energy for the 
intersection of a vortex with a flux tube is [141 
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where m„ is the bare proton mass, m* = m p + 5m* ~ 
m p / 2 [15J is its effective mass, x p ~ 0.05 is the proton 
mass fraction, and pu is the mass density in units of 10 14 
g cm -3 . The range of the interaction between a vortex 
and flux tubes is of order the characteristic radius of a 
flux tube, or the London length, given by 
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The pinning force is F p ~ E p /A*, about 1 MeV fin -1 . 
For a single vortex immersed in a tangle of flux tubes, 
the average length between intersections will equal the 

— 1/2 

average distance between flux tubes, Z$ = , where 

n$ = B/&Q is the areal density of flux tubes, B is the 
average magnetic field strength, and $0 = hc/2\e\ is the 
flux quantum. (This argument is unaltered if flux tubes 
cluster around vortices [9j). The flux tubes are very nu- 
merous: 



?i$ ~ 5 x 10 19 B13 cm 



(7) 



where -B13 is the magnetic field in units of 10 13 G. (The 
relevant field is the average internal field of the star, 
which can be a factor ~ 10 larger than the dipolc field). 
The average separation of pinning intersections along a 
vortex is 



~ 1600B~ 3 1/2 fm 



(8) 



The vortices are far less numerous than the flux tubes, 
- 10 5 (^ s /100 rad s" 1 ) cm" 2 . 

The maximum differential velocity v c between the su- 
perfluid and the flux tube array that can be sustained by 
pinning follows by equating the Magnus force on a vortex 
segment of length to the pinning force per vortex-flux 
tube junction, 



p s Kv c l$ = Ep/A*, 



(9) 



where p s ~ p is the mass density of the neutron super- 
fluid. Combining eqs. [5], [6], [8], and [9], gives [14 1 
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As a result of pinning of vortices to flux tubes, the dif- 
ferential velocity v s approaches v c as the charged compo- 
nent is spun down by the electromagnetic torque, but the 
neutron superfluid is not. Before v c is reached, though, 
vortices will move through thermal activation over the 
pinning barriers presented by the flux tubes. 

Thermally- activated vortex creep. A vortex in the 
outer core spends most of its time pinned to flux tubes. 
Occasionally a thermal excitation frees a segment of the 
vortex, which then moves at the velocity given by eq. 
[2]. After moving a distance Z$, on average, the vortex 
will encounter a flux tube segment to which it pins. The 
probability V that a vortex segment is unpinned at any 
instant is 



V = c- A / T , 



(11) 



where A is the activation energy for unpinning (specified 
below) and T is the temperature in units of energy. The 
average vortex creep velocity is 
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The relaxation time. Assuming axisymmetry, the 
equation of motion of the crust plus charges (component 
c) of spin rate Q c and moment of inertia I c , and the su- 
perfluid (component s) of moment of inertia I s , is given 
by 



di s n s = n c: 



(13) 



where N ext is the external electromagnetic torque, which 
acts only on component c. The angular acceleration of 
the superfluid follows from vorticity conservation [H H(| 

fl s {r, t) = ~ (ittsir, t) + r^Cl 8 {r, t)j v ■ f, (14) 

and is determined by the radial component of the vortex 
velocity which vanishes in the limit of zero dissipation 
(0d = 0). In rotational equilibrium, the crust and super- 
fluid are both spinning down at the same rate Oo, so that 
A'oxt = ^^07 where I = I c + I s . The effects of entrain- 
ment in the core, though responsible for pinning, are not 
very important for the rotational dynamics [14j . and so 
have been ignored here. 

Assuming uniform superfluid rotation dQ s /dr = for 
simplicity, the equation of motion for the angular velocity 
lag cu = f^ s — r2 c is 
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Let the system be in spin equilibrium at t — 0, with Q s — 
fl c = flo for lag loq. An approximate (though accurate) 
solution can be obtained by first expanding A{uj) about 
ujq, to obtain 
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Suppose a glitch occurs at t = 0, through some instability 
that need not be specified for this perturbative analysis. 
Part of the superfluid interior spins down, delivering an- 
gular momentum to the crust and spinning it up; the lag 
lo decreases everywhere. Let the glitch decrease to from 
luq to u)(0+) < ujq. The solution to eqs. [15]-[T7] is 
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where 
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Here r is the intrinsic relaxation time of the system, and 
to is the offset time. The response of the system is de- 
layed by to before relaxing exponentially over a timescale 
r. This behavior is a consequence of the non-linear de- 
pendence of the creep velocity on uj [H, [T(| . 

The activation energy A for a parabolic pinning force 

is 
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where Lo cr = v c /r is the local critical angular velocity lag 
for unpinning [18j |. The equilibrium lag cl>o is given by 
the solution to eq. [TB] with fl s — flo. Introducing the 
spin-down time t s d = f2 c /2|f2o|, the relaxation time from 
eqs. H3, US], and [20], is 



2T 



1/3 



ln- 1/3 [2t 



sd u C r sin 29d 



where the logarithm was evaluated with co 



(21) 
r and 

fl c . The relaxation time increases if E p T, or B are 
increased, and decreases very weakly with the dissipation 
angle 9 d . 

Taking p s — 4 x 10 14 g cm -3 and r — 10 6 cm, the 
result in terms of fiducial values is 
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If fl s is largely unperturbed by the glitch except in a 
localized region, the lag immediately after the glitch will 
be w(0+) = ujq — Af2 c , where Afi c is the observed change 
in the spin rate of the crust. The offset time is then 
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days. (23) 



The total response time to the glitch is ~ to + r, about 
20 days for these fiducial values, and much longer for 
older, cooler stars. For example, for t s d = 10 5 yr, fl c =10 
rad s _1 , -B13 = 1, and T = 3 keV, the relaxation time r 
becomes 18 months. The dependence of t on E p is weak, 
and the relaxation time remains long even if E p =100 
MeV is a significant overestimate. In regions of the star 
where the superfluid spins down, the offset time can be 
considerably larger. Eq. [22] suggests that the coupling 
time in magnetars could be a century or more, but since 
the fields in these objects are comparable to the lower 
critical field for superconductivity(i? c i ~ 10 15 G), the 
pinning description used here is likely to break down. 

Discussion. The chief result of this Letter is that pin- 
ning of vortices to flux tubes in the outer core decou- 
ples the neutron superfluid from the rest of the star over 
timescales of weeks or longer. Given the complexities in- 
volved in measuring post-glitch relaxation, detailed com- 
parison of eq. |21| with individual glitches is perhaps not 
yet warranted, but the calculated response times gener- 
ally agree in magnitude with observed post-glitch relax- 
ation timescales. Since the outer core comprises much 
of the moment of inertia of the star, and much more 
than the inner-crust superfluid, I propose that observed 
post-glitch relaxation represents dynamical response of 
the outer core. For older pulsars, which are also cooler, 
the relaxation time is years. Such long relaxation times 
might explain the nearly step-like nature of many pulsar 
glitches, with little subsequent relaxation. The treatment 
here assumes that at least part of the star is in rotational 
equilibrium before the glitch. Observations show that 
most pulsars do not generally relax completely between 
glitches 1], and this fact should be considered in further 
theoretical analysis and modeling of post-glitch response. 

The results here are at odds with those of Sidery and 
Alpar [llj, who found r ~ 10- 8 (T/10 keV) exp(E p /T), 
giving 10~ 3 (T/10 keV) s for their choice of E p = 0.1 
MeV, a very small pinning energy compared to typical 
estimates of 100 MeV [l^, 14 1. For these values, r given 
by eq. [22] is ~ 2 days. The large discrepancy arises 
from an unphysical expression for the creep velocity used 
in Ref. [n|; see 

A statistical study of the average angular momentum 
transfer in glitches in the Vela and other pulsars [2(| 
shows that glitches are driven by some component of the 



star that has a moment of inertia I g that satisfies 



Ig > 0.016/ C 



(24) 
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recently confirmed with more data 2ll423| . The moment 
of inertia of the inner-crust superfluid is at least several 
percent of the star [24| , supporting the idea that glitches 
arise in the inner crust. Recent work [25| . however, shows 
that most of the neutron mass of the dripped neutrons 
of the inner crust will be entrained by the nuclear clus- 
ters through Bragg scattering with the lattice, an effect 
that occurs whether the neutrons are superfluid or not. 
Entrainment of the neutrons requires the free crust neu- 
trons to account for much more than 1.6% of I c if glitches 
arise in the inner crust, as pointed out in Refs. |22j | and 
[25j . Accounting for entrainment 23 1, eq. [23] becomes 



Ig > 0.08 I C 



(25) 



while the moment of inertia in the inner-crust superfluid 

I,r. IS H 



0.04/, 



(26) 



for a NS of mass 1.4M Q and R — 13 km. If I ~ I c , due to 
tight coupling between the crust and the core superfluid, 
these considerations appear to rule out the inner crust 
superfluid as the origin of glitches. Entrainment causes 
most of the neutron fluid to spin down with the crust, 
and the unentrained conduction neutrons cannot accu- 
mulate angular momentum at a sufficient rate to drive 
large glitches. 

As shown here, however, the effective moment of iner- 
tia of the crust I c represents less of the star than previ- 
ously thought, since the outer core is decoupled from the 
crust over observable timescales. Taking, for example, 
I c = 1/2, cqs. [25] and [26] can both be satisfied, and 
an inner-crust origin of glitches remains viable. Glitches 
might still arise in the core, but an inner-crust origin is 
not excluded as argued in Refs. 22 1 and [23|. Calcula- 
tions of the moment of inertia of the type II region of 
the core are needed to resolve this issue. In any case, the 
long coupling times found in this Letter indicate that the 
core is a viable source of angular momentum for glitches. 

At higher densities in core, the proton condensate is 
expected to change from type II to type I, in which case 
the magnetic flux will be confined to mesoscopic regions, 
probably with a complex geometry. Here vortex pinning 
might be negligible for particular magnetic field geome- 
tries [26j |. On the other hand, in a type I superconductor 
with a complicated field structure, the gradients in the 
magnetic forces could still be large, resulting in vortex 
pinning that might be nearly as effective as that in the 
outer core [27j, leaving almost the entire core decoupled 
by pinning of vortices to magnetic structures. If pinning 
is effective throughout the core, the effective moment of 
inertia of the crust plus charges could be as small as 
~ x p I ~ 0.057, and the angular momentum requirements 
of glitches become considerably easier to satisfy. 

More work is needed to understand the extent to which 
the inner core fluid is coupled to the crust. Finally, I note 



that vortex creep in the outer core is unstable [14j] , which 
could affect the coupling time calculated here. 

I am grateful to N. Chamel for very useful discussions 
about entrainment in the crust. This work was sup- 
ported by NSF Award AST-1211391 and NASA Award 
NNX12AF88G. 



[1 

[2 

[3 

[4 

[5 

[6 

[T 
[8 

[9 

[10; 
[11 
[12 
[13 



[14 

[15, 

[ie; 

[17 
[18 



link@physics.montana.edu; 

C. M. Espinoza, A. G. Lyne, B. W. Stappers, and 
M. Kramer, Mon. Not. Roy. Astron. Soc. 414, 1679 
(2011). 

A. B. Migdal, Nucl. Phys. 13, 655 (1959); V. L. Ginzburg 
and D. A. Kirzhnits, Sov. Phys. JETP 20, 1346 (1965). 

D. J. Dean and M. Hjorth- Jensen, Rev. Mod. Phys. 75, 
607 (2003). 

M. A. Alpar, S. A. Langer, and J. A. Sauls, Astrophys. 
J. 282, 533 (1984). 

M. A. Alpar and J. A. Sauls, Astrophys. J. 327, 723 
(1988). 

I. Easson, Astrophys. J. 228, 257 (1979); M. Abney and 
R. I. Epstein, J. Fluid Mech. 312, 327 (1996). 
P. W. Anderson and N. Itoh, Nature 256, 25 (1975). 
M. A. Alpar, P. W. Anderson, D. Pines, and J. Shaham, 
Astrophys. J. 276, 325 (1984). 

A. D. Sedrakian and D. M. Sedrakian, Astrophys. J. 447, 
305 (2005); A. D. Sedrakian, D. M. Sedrakian, J. M. 
Cordes, and Y. Terzian, ibid. 447, 324 (1995). 
G. Baym, C. Pethick, and D. Pines, Nature 224, 673 
(1969). 

T. Sidery and M. A. Alpar, Mon. Not. Roy. Astron. Soc. 
400, 1859 (2009). 

R. I. Epstein and G. Baym, Astrophys. J. 387, 276 
(1992). 

P. B. Jones, Mon. Not. Roy. Astron. Soc. 253, 279 (1991); 
G. Mendell, Astrophys. J. 380, 515 (1991); H. F. Chau, 
K. S. Cheng, and K. Y. Ding, ibid. 399, 213 (1992); 
M. Ruderman, T. Zhu, and K. Chen, ibid. 492, 267 
(1998). 

B. Link, Mon. Not. Roy. Astron. Soc. 421, 2682 (2012). 
O. Sjoberg, Nucl. Phys. A 65, 511 (1976); N. Chamel 
and P. Haensel, Phys. Rev. C 73, 045802 (2006). 
B. Link, R. I. Epstein, and G. Baym, Astrophys. J. 403, 
285 (1993). 

B. Link and R. I. Epstein, Astrophys. J. 373, 592 (1991). 
This form for the activation energy applies when the vor- 
tex self-energy T v satisfies T v Al / E p l$ << 1. This ratio 
is ~ 0.1 for the outer core. 
[19] The creep velocity used in Ref. [ll[ includes a contri- 
bution for vortex motion against the superfluid flow. To 
move against the flow to a new pinning site a distance 
I® away, a vortex segment of length Z$ must increase its 
energy by an additional amount AE ~ p a K,v c l% sinOd ~ 4 
MeV, for the fiducial values of this Letter. Backward 
transitions are suppressed by a factor ~ e ~ AB / T = e -400 
for T = 10 keV, and so are negligible. As a result, the 
analysis of Ref. [ll[ underestimates the relaxation time 
by a factor ~ 10 12 e~ Ep ^ T . Ref. [ll|] suggested that parts 
of the core could have long relaxation times in their "non- 
linear creep" limit, but did not provide a calculation. 



5 



[20] B. Link, R. I. Epstein, and J. M. Lattimer, Phys. Rev. 

Lett. 83, 3362 (1999). 
[21] R. Dib, V. M. Kaspi, and F. P. Gavriil, Astrophys. J. 

673, 1044 (2008). 
[22] N. Andersson, K. Glampedakis, W. C. G. Ho, and C. M. 

Espinoza, arXiv: 1207.0633 (2012). 



[23] N. Chamel, arXiv: 1210.8177 (2012). 
[24] D. G. Ravenhall and C. J. Pethick, Astrophys. J. 424, 
846 (1994). 

[25] N. Chamel, Nucl. Phys. A 747, 109 (2005); 85, 035801 
(2012). 

[26] A. Sedrakian, Phys. Rev. D 71, 083003 (2005). 
[27] M. Ruderman, Private communication. 



